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1.0  INTRODUCTION 


Radar  backscatter  from  the  ocean  has  been  the  subject  of  considerable 
research  since  World  War  II.  Most  early  studies  aimed  either  at  understanding 
the  nature  of  sea  "clutter"  or  determining  the  clutter  statistics  needed  to 
optimize  radar  parameters  necessary  for  detecting  ships  and  other  targets  in 
the  presence  of  clutter.  These  experimental  programs  indicated  that  radar 
return  from  the  sea  increased  with  wind  speed,  at  least  at  low  wind  speeds, 
but  at  least  one  reported  saturation  at  high  wind  speeds  [Daley,  1973].  This, 
saturation,  however,  was  shown  by  other  investigators  [e.g.,  Claassen,  1972] 
to  be  based  on  misreading  of  the  experimental  results. 

During  the  mid-1960s,  NASA  Johnson  Space  Center  initiated  an  experimental 
program  to  study  the  potential  application  of  radars  to  oceanography.  The 
results  showed  that  a  radar  scatterometer  can  be  used  to  measure  winds  over 
the  ocean.  Based  on  these  results,  Moore  and  Pierson  [1971]  advanced  a 
proposal  for  measuring  world-wide  ocean  winds  with  a  spaceborne  scatterometer. 

Subsequent  measurements  conducted  by  NASA  Langley  Research  Center  under 
the  Advanced  Application  Flight  Experiments  (AAFE)  program  further  verified 
the  dependence  of  scattering  coefficient  on  wind  speed.  A  major  achievement 
of  this  program  was  the  development  by  Jones  et  al .  [1977]  of  circle-flight 
experiments  that  allowed  more  precise  study  of  the  directional  dependence  of 
scattering  from  the  ocean. 

During  the  same  period  SKYLAB  was  launched  and  many  radar  backscatter 
measurements  were  made  with  S-193  scatterometer  from  SKYLAB.  These  measure¬ 
ments  further  confirmed  the  usefulness  of  spaceborne  scatterometry  for  oceanic 
wind  studies. 

SEASAT,  the  first  satellite  dedicated  to  remote  sensing  of  oceans,  was 
launched  by  NASA  in  June  1978.  It  carried  a  number  of  microwave  remote  sen¬ 
sors,  one  of  which  was  a  scatterometer  designed  for  wind  vector  measure¬ 
ments.  A  large  volume  of  data  was  collected  during  the  three  months  of 
satellite  operation,  and  analysis  of  this  data  confirmed  that  a  spaceborne 
scatterometer  can  determine  wind  speed  to  the  required  accuracy. 

The  European  Remote  Sensing  satellite  (ERS-1)  is  scheduled  to  be  launched 
during  1989.  Instrumentation  currently  planned  for  ERS-1  includes  an  active 
microwave  instrument  at  5.3  GHz  designed  to  function  in  three  different 
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modes:  (1)  imaging,  (2)  wave,  and  (3)  wind.  The  dependence  of  radar  cross- 
section  on  wind  speed  at  frequencies  in  C-band  has  not  been  investigated  as 
extensively  as  that  at  higher  frequencies  in  X-  and  Ku-bands.  The  European 
Space  Agency  initiated  an  experimental  program  during  January  -  February  1984 
to  investigate  the  sensitivity  of  radar  backscatter  to  windspeed  at  C-band. 
Tower-based  radar  backscatter  measurements  were  made  by  The  University  of 
Kansas  as  a  part  of  this  experimental  program. 

— rThe  purpose  of  this  report  is  to  summarize  the  data  collected  during 
January  1984  from  the  Nordsee  Tower  in  Germany.  Brief  descriptions  of  the 
system  used  and  the  experiment  are  given  in  Section  2.0.  Angular  responses  of 
the  scattering  coefficient  are  discussed  in  Section  3.0.  The  apparent  ripple 
spectra  calculated  from  backscatter  measurements  are  presented  in  Section  4.0, 
and  the  modulation  transfer  functions  (MTFs)  calculated  from  this  data  set  are 
given  in  Section  5.0.  Section  6.0  presents  some  general  conclusions. 

2.0  SYSTEM  AND  EXPERIMENT  DESCRIPTION 

A  modified  version  of  the  University  of  Kansas  helicopter-borne  spec¬ 
trometer  (HELOSCAT)  was  used  to  acquire  backscatter  data  during  this  experi¬ 
ment.  The  system  measured  radar  scattering  from  the  sea  at  frequencies  from  4 
to  17  GHz  and  incidence  angles  (relative  to  vertical)  from  5°  to  70°. 

'Internal  (relative)  calibration  of  the  system  was  accomplished  by 
measuring  the  signal  from  a  delay  line  switched  in  place  of  the  antenna  before 
and  after  each  data  run.  External  (absolute)  calibration  was  performed  by 
measuring  the  power  received  from  a  target  of  known  radar  cross-section  (a 
metal  sphere  with  a  diameter  of  20  cm).  Normally  a  Luneberg  lens  (which  has  a 
much  higher  cross-section  than  the  metal  sphere)  is  used  for  external  cali¬ 
bration,  but  it  was  not  possible  to  use  the  Luneberg  lens  during  this  experi¬ 
ment  because  of  logistics  problems. x  The  important  specifications  of  the 
system  are  given  in  Table  1.  A  mord  detailed  description  can  be  found  in 
Gogineni  et  al .  [1984]. 

the  main  emphasis  of  these  measurements  was  on  long  data  sets  at  a  few 
fixed  azimuth  angles  (generally  upwind,  downwind  or  crosswind)  and  incidence 
angles  with  VV-polari zation  at  5.3,  10  and  15  GHz.  During  aircraft 
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TABLE  1 

SYSTEM  SPECIFICATIONS 


Type  FM-CW 

Frequency  Range  4-17  GHz 

FM  Sweep  750  MHz 

Transmitted  Power  5-16  dBm 

Intermediate  Frequency  50  kHz 

IF  Bandwidth  13.5  kHz 

Antenna:  Log-Periodic  Reflector 

Polarization  VV,  HH 

Size  61  cm 

Two-Way  Effective  Beamwidths  5.2°,  2.8°  and  2.0°  at 

5.3,  10,  and  15  GHz 

Incident  Angles  18°  to  75°  from  vertical 

Calibration: 

Internal  Delay  Line 

External  Metal  Sphere 

Height  24  meters 


overflights,  measurements  were  made  at  5.3  GHz  at  a  fixed  azimuth  angle  and 
the  same  incidence  angle  used  by  the  systems  in  the  aircraft. 

In  addition,  the  frequency  response  of  ocean  backscatter  was  measured  by 
stepping  the  frequency  at  an  interval  of  400  MHz  in  C-band  and  1  GHz  between  9 
and  17  GHz.  /r  Measurements  at  5.3  GHz  with  VV-polarization  at  a  fixed  azimuth 
angle  were  aTso  made  when  data  with  a  cl  overleaf  buoy  were  being  collected. 

The  experiment  was  conducted  during  January  21  -  February  3,  1984.  The 
windspeed  varied  from  approximately  6  m/sec  to  22  m/sec  during  this  period. 


3.0  RESULTS 

Backscatter  data  from  selected  runs  were  used  to  compute  the  average 
scattering  coefficient  as  a  function  of  wind  speed,  direction  and  incidence 
angle. 
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3.1  Angular  Response 


The  angular  response  of  the  scattering  coefficient  at  5.3  GHz  for 
different  wind  speeds  is  shown  in  Figure  3.1.  The  backscatter  increased  by 
about  1  dB  when  the  wind  speed  increased  from  approximately  7  m/sec  to  13 
m/sec  at  18°.  The  radar  cross-section  difference  between  the  two  wind  speed 
ranges  increased  with  incidence  angle  between  18°  and  40°  and  remained  nearly 
constant  at  angles  greater  than  40°. 

Backscatter  from  the  ocean  near  vertical  is  believed  to  be  dominated  by 
large-scale  roughness.  It  should,  therefore,  decrease  with  increasing  wind 
speed.  The  increase  in  the  scattering  cross-section  with  wind  speed  between 
18°  and  70°  indicates  that  Bragg  scattering  is  the  dominant  mechanism  for 
these  angles  of  incidence. 

Angular  response  of  the  scattering  coefficient  at  5.3,  10  and  15  GHz  is 
shown  in  Figure  3.2.  Backscatter  data  at  5.3  GHz  were  acquired  during 
aircraft  overflights,  and  data  at  10  and  15  GHz  were  collected  shortly  after 
the  overflights  on  January  28,  1984.  At  large  incidence  angles,  the  scat¬ 
tering  coefficient  at  10  GHz  is  about  3-4  dB  higher  than  that  at  5.3  GHz.  The 
increase  in  the  scattering  cross-section  is  less  than  2  dB  between  10  and  15 
GHz.  The  scattering  coefficient  increases  gently  with  frequency. 

The  angular  response  of  the  scattering  coefficient  for  the  upwind 
direction  at  5.3  GHz  is  shown  in  Figure  3.3.  These  measurements  were  made 
during  aircraft  overflights  on  February  3,  1984. 

3.2  Wind  Speed 

One  of  the  most  widely  used  models  to  describe  the  dependence  of  radar 
backscatter  on  wind  speed  is  the  power-law  relationship  and  is  given  by: 

o°(9)  *  a(e)  uy(9) 


where: 

a(0)  =  leveling  coefficient 
y(e)  =  wind  speed  exponent 
u  =  wind  speed  in  m/sec. 
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Scattering  Coofficiant 


FIGURE  3.1: 


Angular  Response  of  Scattering  Coefficient  at  5.3  GHz, 
Crosswind,  VV-Pol . 
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FIGURE  3.2: 


Angular  Response  of  Scattering  Coefficient  at  U  *  6.5  - 
m/sec.  Crosswind,  VV-Pol  .,  January  28,  1984. 
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4.3:  Calculated  Sea  Spectrum  for  Upwind  at  Wind 

15  to  18  m/sec. 


«an 


IQ 


FIGURE  3.10:  Relationship  Between  o°  and  Wind  Speed  at  5.3  GHz 


Scattering  coefficient  vs  windspeed  for  the  upwind  direction  at  5.3  GHz 
and  incidence  angle  of  65°  is  shown  in  Figure  3.10.  The  exponent  at  65°  is 
lower  than  that  at  55°. 


4.0  SEA  SPECTRUM 
4.1  Results 

Selected  radar  backscattering  measurements  at  frequencies  ranging  from 
4.5  -  17  GHz,  and  incidence  angles  ranging  from  45°  to  65°,  were  used  to 
calculate  apparent  ripple  spectra  at  various  windspeeds.  Only  VV  polarization 
was  used  for  these  measurements,  all  of  which  were  in  the  upwind  direction. 

The  results  are  shown  in  Figures  4.1  -  4.4.  In  each  of  these  figures  the  line 
shown  is  a  regression  fit  of  the  form 

W(K)  =  AK"a 

for  the  points  shown.  The  exponent  a  is  tabulated  in  Table  3,  which  also 
includes  values  of  the  coefficient  A. 

TABLE  3 

TABLE  OF  THE  EXPONENTS 


Windspeed  m/s 

a 

A 

9-10 

3.73 

4.07 

X 

10-3 

12  -  13 

3.30 

8.32 

X 

10“3 

15  -  18 

3.03 

7.41 

X 

10'3 

20  -  22 

3.37 

1.380 

X 

10"2 

4.2  Calculation  of  the  Sea  Spectrum,  W(K) 

The  calculation  of  sea  spectrum  values  was  made  using  the  small 
perturbation  theory  [Fung  and  Lee,  1982;  Fung  and  Chan,  1977;  Lawner  and 
Moore,  1985].  According  to  this,  the  backscattering  coefficient  for  vertical 
polarization  is  given  by 


FIGURE  3.9:  Relationship  Between  o°  and  Wind  Speed  at  5.3  GHz 
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FIGURE  3.8:  Relationship  Between  o°  and  Wind  Speed  at  10  GHz 
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FIGURE  3.6:  Relationship  Between  o°  and  Wind  Speed  at  15  GHz 


FIGURE  3.4:  Relationship  Between  o°  and  Wind  Speed  at  5.3  GHz 
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The  relationship  is  generally  expressed  in  log-log  form  as: 


o°(dB)  =  A  +  B  log(u) 


where: 

A  =  10  log  (a(e)) 

B  =  10y 

Scattering  coefficient  as  a  function  of  wind  speed  for  upwind  direction 
at  5.3,  10  and  15  GHz  is  shown  in  Figures  3.4  -  3.6  for  W  polarization.  The 
sensitivity  to  wind  speed  is  higher  at  15  GHz,  but  the  scattering  coefficient 
at  5.3  GHz  also  increases  with  wind  speed.  The  exponent  at  15  GHz  is  much 
higher  than  those  at  5.3  and  10  GHz.  It  is  also  much  higher  than  that  repor¬ 
ted  during  earlier  investigations.  Each  new  experiment  seems  to  produce  expo¬ 
nents  slightly  higher  than  those  during  the  previous  experiments!  Possible 
causes  could  be  either  an  inaccurate  calibration  or  a  loss  of  precision  in 
estimating  the  scattering  coefficient  at  higher  frequencies  because  of  lower 
signal-to-noise  ratio  (at  lower  wind  speeds),  but  the  effect  may  also  be  real. 

The  wind  speed  response  at  5.3  and  10  GHz  for  the  upwind  direction  with 
HH-polarization  is  shown  in  Figures  3.7  and  3.8,  respectively.  The  sensitiv¬ 
ity  to  wind  speed  at  5.3  GHz  is  higher  than  that  at  10  GHz.  The  data  reported 
in  these  figures  were  collected  during  different  times.  The  scatter  in  the 
data  points  at  10  GHz  might  be  because  of  a  change  of  environmental  parameters. 

Wind  speed  response  at  5.3  GHz  for  the  crosswind  direction  is  shown  in 
Figure  3.9.  The  exponent  (2.15)  for  crosswind  direction  is  lower  than  that 
for  upwind  (2.67)  direction.  (See  Table  2). 


TABLE  2 


WINDSPEED  EXPONENTS 


Figure 

Freq. 

(GHz) 

Pol . 

Wind  Dir. 

Inc.  Angle 
(°) 

104a 

7 

3.4 

5.3 

VV 

upwind 

55 

.156 

2.67 

3.5 

10  VV 

upwind 

55 

.234 

2.31 

3.6 

15  VV 

upwind 

55 

.021 

3.79 

3.7 

5.3 

HH 

upwind 

55 

.215 

2.37 

3.8 

10  HH 

upwind 

55 

1.271 

1.81 

3.9 

5.3 

VV 

crosswind 

55 

.386 

2.15 
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FIGURE  3.3: 


Angular  Response  of  Scattering  Coefficient  at  5.3  GHz. 
Upwind,  VV-Pol .,  February  3,  1984.  Aircraft  Overflights. 


(4.1) 


<JVV(M)  =  8k4  1^  I2  S(K,4>)/K 

where  K  is  the  ocean-surface  wavenumber,  and  $  is  an  azimuthal  coordinate.  We 
defined  the  sea  spectrum,  W(K)  as  follows: 


S(K,») 

K 


(4.2) 


For  the  VV  case,  a  good  approximation  for  the  average  radar  backscattering 
coefficient  is  [Fung  and  Lee,  1983] 


8k  |auur  W(K) 


(4.3) 


where: 


k  =  radar  wavenumber 
K  =  2sin8  =  ocean  wavenumber 
0  =  incident  angle. 

£  —  1 

avv  *  Rvcos2e  +  1-fe - )  U  +  Ry)2  sin2® 


epcos0  -  (£p  -  sin29)^2 

V  epcos0  +  (e~  -  sin20)^2 

ep  =  relative  complex  dielectric  constant. 

Thus  calculation  of  the  sea  spectrum,  W(K),  was  made  using  the  equation 


8k  i°vvr 


(4.4) 


4.3  Comparison  with  Previous  Measurements  and  with  Model 

The  results  obtained  from  this  experiment  can  be  compared  with  the  semi- 
empirical  model  proposed  by  Fung  and  Lee,  [1982]/  For  comparison,  this  model 
predicts  that,  for  the  range  of  K  shown  in  our  figures,  the  value  of  W(K)  is 
given  by  [Lawner  and  Moore,  1985] 


21 


W(K) 


where: 


.875  (2it)p"2(1  +  ^-)  g 

,,2  #p+l» 

K[K(1  +  2  ' 

kl 


km*aT*  (3.73) 

.2 


2  -2 
cm 


g  =  981  cm/s 
p  =  sea  water  density 
t  =  sea  surface  tension 


U  * 2  tw» 


and  R  is  the  ratio  of  slope  variances  in  the  crosswind  and  upwind 
directions  as  given  by  Cox  and  Munk's  [1954]  clean  sea  model: 


ct 


.003  +  1.92  x  10"3  u 


12.5 


3.16  x  10"3  u. 


"ut  . 12.5 

and  u12  5  =  wind  speed  at  a  height  of  12.5  meters  in  meters  per  second. 


The  results  shown  in  Figures  4.1  -  4.4  fall  much  closer  to  the  semi- 
empirical  model  than  earlier  tower  results  made  in  1979  from  the  Noordwijk 
tower  [Lawner  and  Moore,  1985],  The  earlier  results  could  only  be  brought 
into  good  agreement  with  the  model  after  compensating  empirically  for  the 
difference  between  tower-based  and  air/spacecraft  measurements.  This  seemed 
reasonable  since  the  model  was  based  on  air/spacecraft  measurements  as  opposed 
to  tower-based  measurements.  However,  the  present  results  do  not  seem  to  need 
such  an  adjustment  since  they  are  already  in  fairly  close  agreement  with  the 
model.  The  reason  for  this  that  there  is  little  bias  between  the  tower-based 
and  aircraft  measurements  at  large  incidence  angles. 


5.0  MEASUREMENTS  OF  THE  MODULATION  TRANSFER  FUNCTION 

Consider  a  microwave  antenna  illuminating  a  portion  of  the  ocean  surface 
small  in  comparison  with  the  wavelength  of  any  long  ocean  waves.  The  primary 
scatterers  are  short  waves  (ripples),  a  few  centimeters  in  length,  that  are 
riding  on  the  long  ocean  waves.  Since  the  radar  only  detects  the  short  waves. 
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the  detectability  of  the  long  ocean  waves  by  radar  becomes  possible  by  modu- 
latng  the  radar  return  power  by  the  long  ocean  waves. 

For  low  to  moderate  sea  states,  the  linear  model  introduced  by  Keller  and 
Wright  [1975]  is  used  to  describe  the  modulation.  As  shown  in  Figure  5.1,  the 
proposed  model  assumes  that  the  scattered  power  is  linearly  dependent  on  the 
long  ocean  wave  through  a  transfer  function  R(f). 


Using  linear  network  theory,  with  the  height  of  a  long  ocean  wave  represented 
by  W(t)  and  radar  return  power  by  P(t),  we  find 


R(f) 


(5.1) 


where  Sw  w  is  the  spectrum  of  the  waveheight  and  Sy  p  is  the  cross-spectrum  of 
the  wave  height  and  power.  Note  that  the  above  equation  can  be  represented  in 
the  time  domain  by 


P(t)  =  V  +  RW(t) 


where  V  is  the  mean  scattered  power. 

Normalizing  (5.1)  by  l5"  and  defining 


9 


(5.2) 


the  modulation  transfer  function  m(f),  as  defined  by  Keller  and  Wright  [1975] 
is  obtained. 
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m(f)  = 


3W.P 


kTs 


U,W 


(5.3) 


where  K  is  the  wave  number  of  the  long  ocean  waves  at  frequency  f. 

An  Important  quantity  in  analysis  of  this  type  is  the  coherence  function 
defined  by 


coh2(f)  = 


If  Sp>p  is  proportional  to  SWjW,  the  coherence  function  has  a  value  of  one  at 
any  given  frequency.  Values  less  than  one  indicate  either  (a)  noisy  signals, 
(b)  nonlinearities,  or  (c)  Sp  p  is  due  to  sources  in  addition  to  Sw  w. 

Two  sources  of  modulation  of  the  scattered  microwave  power  exist.  These 
are  the  tilt  modulation  and  the  hydrodynamic  modulation.  Therefore,  the 
modulation  transfer  function  can  be  broken  down  into  two  components: 


R  “  Rtilt  +  Rhydro 


(5.4) 


5.1  Tilt  Modulation 

The  tilt  modulation  is  due  to  the  fact  that  the  ripples  are  seen  by  the 
radar  at  different  local  Incidence  angles  depending  on  their  location  on  the 
long  waves.  The  backscatter  from  the  ripples  involves  Bragg  resonance  with  a 
ripple  component  of  wavelength  Ag  =  A/2  s1n0^,  where  A  is  the  radar  wavelength 
and  9L  is  the  local  angle  of  incidence. 

The  Bragg  scattering  cross-section  per  unit  area  for  a  tilted  facet  is 
given  by 


«Sp(0Ll*)  =  8K4a14cos40L  |app|2  W4(0ulO 


(5.5) 


Assuming  that  the  short  wave  spectrum  Is  a  Phillips  spectrum,  (5.5)  takes  the 
form 
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-  [5  (AK)3  ♦ 


16  °5(AK)' 


...J 


sin3<j>  + 


(5.9) 


Notice  that  harmonics  of  the  ocean  wave  frequency  are  present  in  the 
scattered  power  and  the  amplitude  of  each  harmonic  is  nonli nearly  related  to 
the  ocean  waveheight.  A  measure  of  first-order  amplitude  nonlinearity  can  be 
obtained  by  considering  the  ratio 


r 


1 


3/4  a3(AK)3 
(AK)ax 


=  |^(ak)2 


(5.10) 


Similarly,  for  the  lowest-order  harmonic  nonlinearity  we  can  define  r 2  as 
«?(AK)2 

r2 =  -tkrsj — =  "  1/2  (AK)  ( }  (5*U) 


The  ratios  r^  and  r2  at  different  pointing  angles  and  for  slopes  of  0.1 
and  0.3  are  given  in  Table  4.  As  seen,  at  specific  pointing  angles  and  for 
steep  slopes  of  the  long  waves,  the  process  is  nonlinear.  One  consequence  of 
the  nonlinearities  is  that  the  amplitude  of  the  tilt  modulation  transfer 
function  decreases  as  the  slope  increases. 


TABLE  4 

RATIOS  rx  AND  r2  AT  DIFFERENT  SLOPES  AND  POINTING  ANGLES 


9P 

r^AK-O.l) 

r1(AK=0.3) 

r2(AK=0.1) 

r2(AK=0.3) 

20 

0.320 

2.880 

0.371 

1.110 

30 

0.152 

1.371 

0.259 

0.777 

40 

0.096 

0.864 

0.212 

0.636 

50 

0.075 

0.675 

0.194 

0.582 

60 

0.071 

0.639 

0.202 

0.606 

70 

0.100 

0.901 

0.251 

0.753 

To  Illustrate  the  effects  of  the  nonlinearities,  we  assume  that  only  the 
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first  three  terms  of  (5.7)  are  significant.  Then  we  have 


P(t) 


a  +  an  S  +  a,S 
0  1  i 


Assuming  that  slopes  are  normally  distributed,  it  is  shown  in  Appendix  1  that 
the  spectrum  of  the  received  power  is  given  by: 

Spp(w)  «  c ' 6(uj)  +  «12“2SWjW  +  7  (0t10t2  +  a2^  *  (“\w>J 

Also  in  Appendix  1,  the  cross-spectrum  of  wave  height  and  power  has  been 
derived  as 


SW,P(t^  = 


a^u) 


SW,W  +  a2^  ^“2sW,W^  *  ^2sW,W^ 


where  *  denotes  convolution. 


Then,  from  (5.3)  we  have 


m(f) 


a2(l/2.)[(»2St,  w)  ♦  (.2SH_H)j 
KP  «JkRm  ’ 


Notice  that  the  second  term  of  the  above  equation  is  due  to  the  effects  of  the 
nonlinearities. 


5.2  Hydrodynamic  Modulation 

If  one  assumed  that  the  ripple  amplitudes  were  uniform  over  a  large-scale 
wave,  the  modulation  of  the  received  signal  would  be  completely  governed  by 
the  tilt  modulation.  Hydrodynamic  modulation  of  the  ripples  causes  a  non- 
uniform  distribution  of  the  ripples  over  the  long  ocean  waves. 

If  the  spectrum  of  the  short  waves  S(K)  and  the  slopes  of  the  long  waves 
are  known,  the  scattering  coefficient  due  to  the  tilt  effect  can  be  computed 
as  follows: 

3 

°VV(t)  =  iflnY  lRvcos2eL  +  1/2(1  “  e){1  +  Rv)2s1n20L|2s(6ksin6L)  (5*12) 
where: 

®L  *  9P  +  Tan_1S 
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Ry  3  Fresnel  reflection  coefficient 
e  =  dielectric  constant  of  the  sea 


S(*)  »  sea  spectrum 

The  tilt  modulation  can  be  removed  from  the  scattering  coefficient  to  obtain 
Ao°(t)  3  o°  (t)  -  o°h  (t)  (5.13a) 

0  ,  V 

<4  (t)  »  °m  -  ■  4  (t)  (5.13b) 

8VV  (t) 

where  <r  (t)  is  the  measured  backscattering  coefficient,  and  o°v  (t)  is  given 
in  (5.12).  Normalization  shown  in  (5.13b)  is  necessary  because  the  absolute 
level  of  the  theoretical  o°  used  for  tilt  modulation  may  not  agree  with 
experiment.  A  hydrodynamic  modulation  transfer  function  can  then  be  calcula¬ 
ted  by  using  the  modified  scattering  coefficient,  instead  of  the  actual  scat¬ 
tering  coefficient,  in  the  MTF  calculations. 

5.3  Calculation  of  the  MTF 

Equation  (5.1)  is  used  to  calculate  the  MTFs.  The  spectral  quantities  in 
(5.1)  must  be  replaced  by  the  weighted  sum  of  their  periodograms,  obtaining 
the  following  equation: 


.  mponi  ,5.14) 

l«(f)l2 

where  the  overbar  denotes  ensemble  averaging.  Accordingly,  the  coherence 
function  can  be  calculated  as 

coh2(f)  3  -  lP{f)WV) Z-  (5.15) 

|P(f) I2  |W(f) | 2 

In  order  for  the  periodograms  in  the  above  equation  to  be  good  estimates 
of  the  power  spectral  density,  averaging  over  many  subsets  of  a  long  data  set 
is  required.  This,  in  turn,  requires  that  the  environmental  conditions  over 
the  whole  period  of  the  run  be  reasonably  stable.  Unfortunately,  this  has  not 
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been  the  case  for  the  Nordsee  data.  Figure  5.2  shows  the  wave  height  spectra 
of  two  seven-minute  periods  separated  by  a  ten-minute  interval.  Comparing  the 
two  power  spectra,  notice  that  the  waveheight  spectrum  in  the  second  graph  has 
shifted  to  the  right.  This  shift  in  the  dominant  waveheight  frequency  does 
not  allow  averaging  over  the  entire  period  to  obtain  more  degrees  of  freedom. 
Considering  these  problems,  the  long  runs  were  divided  into  subsets  of  100 
seconds  each  to  calculate  the  spectra.  Each  subset  was  padded  with  zeros  to 
make  a  128-second  subset.  The  Fourier  transform  of  each  subset  was  multiplied 
by  a  moving  average  window  function.  After  multiplying  by  the  window  func¬ 
tion,  seven  128-second  runs  were  averaged  to  compute  the  associated  power 
density  spectra. 

With  the  environmental  conditions  changing  rapidly,  it  was  necessary  to 
shorten  the  run  time.  However,  use  of  the  shortened  time  span  to  accommodate 
the  rapidly  changing  environmental  conditions  caused  the  spectra  to  have  large 
vari ances. 

To  increase  the  accuracy  of  the  spectra,  a  more  advanced  method  of  spec¬ 
tral  estimation  is  needed.  The  Singular  Value  Decomposition  (SVD)  method 
introduced  by  Cadzow  [1982]  is  probably  the  most  suitable  method.  A  program 
currently  under  development  at  the  Remote  Sensing  Laboratory  of  the  University 
oi  sansas  employs  this  method  to  compute  the  power  spectral  densities  and 
MTFs.  Upon  completion  of  the  testing  of  the  testing  of  this  program*  it  will 
be  used  in  computing  the  MTFs. 


5.4  Preliminary  Results 

Figure  5.3  shows  a  sample  run  in  the  upwind  direction  at  a  radar  fre¬ 
quency  of  10  GHz  and  a  wind  speed  of  18.8  m/s.  The  spectrum  of  the  wave 
height  shows  a  single  peak  at  a  frequency  of  0.15  Hz. 

The  magnitude  and  phase  of  the  MTF  are  plotted  in  Figures  5.3b  and  5.3c, 
respectively.  The  phase  of  the  MTF  is  close  to  zero  at  the  point  of  high 
coherency  and  the  MTF  has  a  magnitude  of  about  0.8  at  that  frequency.  The  low 
values  of  the  MTF  clearly  illustrate  that  at  high  wind  speeds  (high  slopes) 
the  assumed  linear  relationship  is  weak  and  the  process  could  be  better 
described  by  a  nonlinear  model. 
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FIGURE  5.2:  Spectra  of  Wavehelght  Calculated  from  Vertically  Polarized  5.3 

GHz  Backscatter  at  a  Pointing  Angle  of  55°  and  at  a  Wind  Speed 
of  8.8  n/s  In  the  Upwind  Direction. 
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One  method  of  improving  the  calculation  is  to  calculate  the  MTFs  using 
the  A o°,  as  explained  in  the  previous  sections.  The  result  is  shown  in  Figure 
5.4.  Notice  that  the  magnitude  of  MTF  has  increased.  This  demonstrates  the 
fact  that  removing  the  tilt  modulation  to  some  extent  has  linearized  the  pro¬ 
cess.  The  higher  values  of  coherence  function  and  cross-correlation  in  Figure 
5.4  agree  with  this  conclusion.  Since  the  MTF  in  Figure  5.4  is  the  hydrody¬ 
namic  portion  of  the  modulation  transfer  function,  one  could  argue  that  the 
hydrodynamic  modulation  is  more  linear  than  the  tilt  modulation  at  high  wind 
speeds. 

Figure  5.5  shows  the  spectrum  of  the  waveheight,  MTF,  phase  of  the  MTF, 
coherence  function  and  the  cross-correlation  function  for  vertically  polarized 
15  GHz  backscatter  at  a  wind  speed  of  16.5  m/s  in  the  upwind  direction.  The 
magnitude  of  the  MTF  is  approximately  equal  to  the  magnitude  of  the  MTF  for 
the  10  GHz  case,  but  the  phase  is  close  to  320°  in  this  case. 

Figure  5.6  shows  the  plots  of  hydrodynamic  modulation  transfer  function 
and  the  cross-correlation  of  Ac0  and  waveheight.  Notice  that  the  MTF  has 
decreased  compared  to  the  previous  case  while  its  phase  has  essentially 
remained  unchanged. 

In  conclusion,  we  emphasize  that,  since  the  microwave  power  is  nonlin- 
early  related  to  the  ocean  wave  height,  the  concept  of  the  MTF  is  not  appli¬ 
cable  for  higher  waves  and  slopes.  A  nonlinear  model  is  needed  to  describe 
the  relationship  between  the  scattered  power  and  the  ocean  wave  height. 


6.0  CONCLUSIONS 

This  paper  extends  the  results  presented  earlier  to  include  the  sea 
spectrum  computed  from  the  scattering  coefficients  and  the  modulation  transfer 
function  computed  from  the  selected  data.  In  addition,  problems  associated 
with  the  computation  of  the  modulation  transfer  functions  are  discussed.  The 
scattering  coefficients  presented  are  accurate  up  to  +  1  dB  and  no  systematic 
calibration  errors  have  been  observed. 

It  was  not  possible  to  make  comparisons  with  the  aircraft  data  because 
they  are  not  presently  available  to  The  University  of  Kansas.  A  detailed  data 
analysis  will  have  to  await  receipt  of  results  from  the  other  investigators. 
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FIGURE  5.4:  Sample  Run  for  a  Vertically  Polarized  10  GHz  Backscatter  at  a 

Pointing  Angle  of  64.8°  and  at  a  Hind  Speed  of  18.8  m/s  in  the 
Upwind  Direction. 


FIGURE  5.5:  Sample  Run  for  a  Vertically  Polarized  15  GHz  Backscatter  at  a 

Pointing  Angle  of  55°  and  at  a  Hind  Speed  of  16.5  m/s  In  the 
Upwind  Direction. 


FIGURE  5.6:  Sample  Run  for  a  Vertically  Polarized  15  GHz  Backscatter  at  a 

Pointing  Angle  of  55°  and  at  a  Wind  Speed  of  16.5  m/s  In  the 
Upwind  Direction. 
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APPENDIX  1 

DERIVATION  OF  THE  RADAR  SIGNAL  SPECTRA 
Assume  that  P(t)  *  aQ  +  S  +  S2.  Then 

Rpp(T)  a  E  {[°0  +  OjSU)  +  a2  S2(t)j[aQ  +  S(t+x)  +  ^  S2)t+t)J} 

=  aQ2  +  E  (S(t)S(t+T)}  o^2  +  E[S(t+r)J 

+  V2  £[s2(t+T)J  +  aoal  EtS(tH  +  ala2  S2(t+T)j 

+  aQa2  E  [  S2  ( t )  J  +  E[S2(t)  S(t+x)j  +  a*  E[S2(t)S2(t+r)  J 

Assume  that  the  slopes  are  normally  distributed  with  mean  equal  to  zero.  Let 
S(t)  =  S  and  Y  =  $(t+T) 

Then 

E  [ S ( t )  S2(t+r) J  *  E  [SY2J 

r2o2 

E  [SY2 j  =  E  [S*E{Y2|S}J  =  E  {S[a2  (1-r2)  +— pX2J} 

al 

2  2  2  2 

where  =  E(S  )  and  o2  =  E(Y  )  and  r  is  the  correlation  coefficient.  Thus, 

E(SY2)  E(S3)  .  .  E2CYS] 

°1  ®1 

Therefore, 

E  [ S ( t )  S2(t+x)J  =  E2[S(t)  S(t+x)j  =  R*$  (T) 

The  autocorrelation  of  the  power  can  be  simplified  to 
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